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Due to the ultrarelativistic velocity of the expanding "fireshell" (Lorentz gamma factor 
7 ~ 10^ — 10'^), photons emitted at the same time from the fireshell surface do not 
reach the observer at the same arrival time. In interpreting Gamma-Ray Bursts (GRBs) 
it is crucial to determine the properties of the EQuiTemporal Surfaces (EQTSs): the 
locus of points which are source of radiation reaching the observer at the same arrival 
time. In the current literature this analysis is performed only in the latest phases of the 
afterglow. Here we study the distribution of the GRB bolometric luminosity over the 
EQTSs, with special attention to the prompt emission phase. We analyze as well the 
temporal evolution of the EQTS apparent size in the sky. We use the analytic solutions 
of the equations of motion of the fireshell and the corresponding analytic expressions of 
the EQTSs which have been presented in recent works and which are valid for both the 
fully radiative and the adiabatic dynamics. We find the novel result that at the beginning 
of the prompt emission the most luminous regions of the EQTSs are the ones closest to 
the line of sight. On the contrary, in the late prompt emission and in the early afterglow 
phases the most luminous EQTS regions are the ones closest to the boundary of the 
visible region. This transition in the emitting region may lead to specific observational 
signatures, i.e. an anomalous spectral evolution, in the rising part or at the peak of the 
prompt emission. We find as well an expression for the apparent radius of the EQTS in 
the sky, valid in both the fully radiative and the adiabatic regimes. Such considerations 
are essential for the theoretical interpretation of the prompt emission phase of GRBs. 

Keywords: gamma rays: bursts — relativity 

1. Introduction 

It is widely accepted that Gamma-Ray Burst (GRB) afterglows originate from the 
interaction of an ultrarelativistically expanding shell into the CircumBurst Medium 
(CBM). Differences exists on the detailed kinematics and dynamics of such a shell 
(see e.g. Refs. [U [2]and Refs. therein). 

Due to the ultrarelativistic velocity of the expanding shell (Lorentz gamma factor 
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7 ~ 10^ — 10^), photons emitted at the same tmie in the laboratory frame (i.e. the 
one in which the center of the expanding shell is at rest) from the shell surface but 
at different angles from the line of sight do not reach the observer at the same arrival 
time. Therefore, if we were able to resolve spatially the GRB afterglows, we would 
not see the spherical surface of the shell. We would see instead the projection on the 
celestial sphere of the EQuiTemporal Surface (EQTS), defined as the surface locus 
of points which are source of radiation reaching the observer at the same arrival 
time (see e.g. Refs. ElllElEllIllHllllllOland Refs. therein). The knowledge of the 
exact shape of the EQTSs is crucial, since any theoretical model must perform an 
integration over the EQTSs to compute any prediction for the observed quantities 
(see e.g. Refs. |lllll2l[9l[inill3l[2l[Illll5]and Refs. therein). 

One of the key problems is the determination of the angular size of the visible 
region of each EQTS, as well as the distribution of the luminosity over such a visible 
region. In the current literature it has been shown that in the latest afterglow phases 
the luminosity is maximum at the boundaries of the visible region and that the 
EQTS must then appear as expanding luminous "rings" (see e.g. Refs. [HI [5J [SJ [71 
[TTlfTBlfroil^l^l^ and Refs. therein). Such an analysis is a pplied only the late st 
afterglow phases to interp ret data from radio observations l^'j|to|lij|zl|id | ^4|^o | 
gravitational microlensing | 26 | 27 | 28 [29j rpj^^ shell dynamics is usually assumed to 
be fully adiabatic and to be described by a power-law 7 cx r~^/^, following the 
Blandford-McKee self similar solution'^, where 7 and r are respectively the Lorentz 
gamma factor and the radius of the expanding shell. Such a power-law behavior has 
been extrapolated backward from the latest phases of the afterglow all the way to 
the prompt emission phase. 

In Refs. l9l [TOl [T] there have been presented the analytic solutions of the equations 
of motion for GRB afterglow, compared with the above mentioned approximate 
solutions, both in the fully radiative and adiabatic regimes, and the corresponding 
analytic expressions for the EQTSs. It has been shown that the approximate power- 
law regime can be asymptotically reached by the Lorentz gamma factor only in the 
latest afterglow phases, when 7^-^10, and only if the initial Lorentz gamma factor 70 
of the shell satisfies 70 > lO'^ in the adiabatic case or 70 > 10"* in the radiative case. 
Therefore, in no way the approximate power-law solution can be used to describe the 
previous dynamical phases of the shell, which are the relevant ones for the prompt 
emission and for the early afterglow. 

Starting from these premises, in this Paper we present the distribution of the 
extended afterglow luminosity over the visible region of a single EQTSs within the 
"fireshell" model for GRBs. Such a model uses the exact solutions of the fireshell 
equations of motion and assumes a fully radiative dynamics (see Refs. [3TJ [32| and 
Refs. therein for details). We recall that within the fireshell model the peak of the 
extended afterglow encompasses the prompt emission. We focus our analysis on the 
prompt emission and the early afterglow phases. Our approach is therefore comple- 
mentary to the other ones in the current literature, which analyze only the latest 
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afterglow phases, and it clearly leads to new results when applied to th e pr ompt 

lOOl 

emission phase. For simplicity, we consider only the bolometric luminosity^, since 
during the prompt phase this is a good approximation of the one observed e.g. by 
BAT or GBM instrumentJSSMI. 

The analysis is separately performed over differ- 
ent selected EQTSs. The temporal evolution of the luminosity distribution over the 
EQTSs' visible region is presented. As a consequence of these results, we show the 
novel feature that at the beginning of the prompt emission the most luminous re- 
gions of the EQTSs are the ones closest to the line of sight. On the contrary, in the 
late prompt emission and in the early afterglow phases the most luminous EQTS 
regions are the ones closest to the boundary of the visible region. This transition in 
the emitting region may lead to specific observational signatures, i.e. an anomalous 
spectral evolution, in the rising part or at the peak of the prompt emission. We 
also present an analytic expression for the temporal evolution, measured in arrival 
time, of the apparent radius of the EQTSs in the sky. We finally discuss analogies 
and differences with other approaches in the current literature which assumes an 
adiabatic dynamics instead of a fully radiative one. 



2. The Equitemporal Surfaces (EQTS) 

For the case of a spherically symmetric fireshell considered in this Letter, the EQTSs 
are surfaces of revolution about the line of sight. The general expression for their 
profile, in the form "d = "d{r), corresponding to an arrival time ta of the photons at 
the detector, can be obtained from (see e.g. Ref. fTO)) : 

eta = ct (r) — r cos -f r* , (1) 

where r* is the initial size of the expanding fireshell, d is the angle between the 
radial expansion velocity of a point on its surface and the line of sight, t = t{r) is its 
equation of motion, expressed in the laboratory frame, and c is the speed of light. 

In the case of a fully radiative regime, the dynamics of the system is given by 
the following solution of the equations of motion (see e.g. Refs. |35l [10] and Refs. 
therein): 

^ 1 + (Mcbm/A^B) (1 + 7o~') [1 + (1/2) (Mebm/A/j?)] 
^ lo^ + {M.^njMB) (1 + 7o"') [1 + (1/2) (Mcbm/Ms)] ' 

where 7 is the Lorentz gamma factor of the fireshell, Mcbm is the amount of CBM 
mass swept up within the radius r and 70 and Mb are respectively the values 
of the Lorentz gamma factor and of the mass of the fireshell at the beginning of 
the extended afterglow phase. Correspondingly, the exact analytic expression for 
t = t{r) iPl 



(3) 

/C f [A+ir/r^)f{A-^ + l) \ m>„ fr^-rj) 

12cm°A^ ]^ [A3 + (r/ro)3](A+l)3 J + ^° + ScVC J 
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where A = {/{Mb - m°)/m°, C = Ms'^ilo - l)/(7o + 1), to is the value of the 
time t at the beginning of the extended afterglow phase and m° — (4/3)7rmpncbm''o- 
Inserting Eq.® into Eg. ([T]) we have the analytic expression for the EQTS in the 
fully radiative regimel^ 



„o«,9 - (Mfl-O('-'-o) 



arctan 



2r—Aro 



arctan 



2-A 
A^/3 



o / 4 4 \ 

SrVC V '1 J 



In 

12rm?A2 1 \A^ + (r/r^f\{A+lf 
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Instead, the corresponding equations in the adiabatic regime ar^^^ 

2 _ 7o + (M.^JMb) + (Afebm/Ms)' 



cos 



1 + (Afebm/Ms) + (Afcbm/Az/ij)' 



(4) 



(5) 

(6) 

(7) 




Fireshell radius (cm) 

Fig. 1. The Fireshell Lorentz gamma factors in the fully radiative regime (red line), given by 
Eq-©, and in the adiabatic regime (blue line), given by Eq. Jsjl. 

A comparison between the Lorentz gamma factors in the two regimes is presented 
in Fig. [H Here and in the following we assume the same initial conditions as in 
Ref.Iini namely 7o = 310.131, ro = 1.943 xlO^'' cm, =6.481x10^8, r* = 2.354 x 
10^ cm, ricbm = 1.0 particles/cm'^, Mb = 1.61 x lO'^" g. For simplicity, and since we 
are interested in the overall behavior of the luminosity distribution, we assume a 
constant CBM density, neglecting the inhomogeneities which are responsible of the 
temporal variability of the prompt emission!^. 
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3. The extended afterglow luminosity distribution over the EQTS 
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Fig. 2. We plot the luminosity density C given in Eq-JS} over 8 different EQTSs, corresponding to 
8 different ta values ranging from the prompt emission (5 seconds) to the early (1 hour) afterglow 
phases (see title above each plot). Each EQTS is represented as its projection on a plane containing 
the line of sight, which coincides with the X axis of each plot and which is represented by a dashed 
gray line. The observer is far away along the X axis. The luminosity density distribution over each 
EQTS is represented by a color gradient, with the highest values corresponding to the lightest 
colors (see the key on the left of each plot). The scales of the different axes and of the color 
gradient are different among the different plots, since it was not possible to choose a single scale 
suitable for all of them. The black curves represent the condition cos ^ = v/c, see Eq.®. The 
boundaries of the visible regions are therefore defined by the intersections of such lines with the 
EQTS external profiles, which coincide with the points where the EQTS external profiles show an 
horizontal tangent in the plots. 
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Within the firesheh model, the GRB extended afterglow bolometric luminosity 
in an arrival time dta and per unit solid angle dVl is given by (details in Refs. 1331 
l32land Refs. therein): 



where Ae is the energy density released in the interaction of the ABM pulse with the 
CBM measured in the comoving frame, A = 7(1 — {v/c) cosd) is the Doppler factor, 
dYi is the surface element of the EQTS at arrival time ta on which the integration 
is performed, and it has been assumed the fully radiative condition. We are here 
not considering the cosmological redshift of the source, which is constant during the 
GRB explosion and therefore it cannot affect the results of the present analysis. We 
recall that in our case such a bolometric luminosity is a good approximation of the 
one observed in t he prom pt emission and in the early afterglow by e.g. the BAT or 
GBM instrumentPSEl. 

We are now going to show how this luminosity is distributed over the EQTSs, i.e. 
we are going to plot over selected EQTSs the luminosity density £(r, ip; ta). The 
results are represented in Fig. [2] We chose eight different EQTSs, corresponding 
to arrival time values ranging from the prompt emission (5 seconds) to the early 
(1 hour) afterglow phases. For each EQTS we represent also the boundaries of 
the visible region due to relativistic collimation, defined by the condition (see e.g. 
Refs. l36l and Refs. therein): 



We obtain that, at the beginning (ta — 5 seconds), when 7 is approximately con- 
stant, the most luminous regions of the EQTS are the ones along the line of sight. 
However, as 7 starts to drop {ta ^ 30 seconds), the most luminous regions of the 
EQTSs become the ones closest to the boundary of the visible region. This tran- 
sition in the emitting region may lead to specific observational signatures, i.e. an 
anomalous spectral evolution, in the rising part or at the peak of the prompt emis- 
sion. 

4. The EQTS apparent radius in the sky 

From sec. [3] we obtain that within the fircshcll model the EQTSs of GRB extended 
afterglows should appear in the sky as point-like sources at the beginning of the 
prompt emission but evolving after a few seconds into expanding luminous "rings" , 
with an apparent radius evolving in time and always equal to the maximum trans- 
verse EQTS visible radius r± which can be obtained from Eas. (|ll9p : 





cosi^ > v/c. 



(9) 




(10) 
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where t = t{r) is given by Eq.(|3]). With a small algebra we get: 

1 = rhir) 



ta = t[r)-{T/c)^l--i{r)-^ + {r*/c) ' ^^^^ 

where 7 = 7 (r) is given by Eq.Q and t = t(r) is given by Eq.Q, since we 
assumed the fully radiative condition. Eq. (fTT|) defines parametrically the evolution 
of r_L = r_L {ta)^ i.e. the evolution of the maximum transverse EQTS visible radius 
as a function of the arrival time. In Fig. [2] we saw that such rx coincides with the 
actual value of the EQTS apparent radius in the sky only for ta ^ 30 s, since for 
^ 30 s the most luminous EQTS regions are the ones closest to the line of sight 
(see the first three plots in Fig. [2]). Therefore, in Fig. |3]we plot r±_ given by Eq. (fTT|) 
in the fully radiative regime together with the actual values of the EQTS apparent 
radius in the sky taken from Fig. [2] in the three cases in which they are different. 
It is clear that, during the early phases of the prompt emission, even the "exact 
solution" given by Eq.(fTT|) can be considered only an upper limit to the actual 
EQTS apparent radius in the sky. 




Max. trans, visible radius, fully radiative 

Actual value 



10^ 10" 
Arrival time (tg) (s) 



10^ 



Fig. 3. The EQTS apparent radius in the sky as a function of the arrival time ta in the fully 
radiative regime. The red line represents the maximum transverse EQTS visible radius following 
Eg. mil l. The three blue points represents the actual EQTS apparent radius in the sky in the three 
cases ta = 5.0 s, ta = 10.0 s and ta = 30.0 s in which, as shown in Fig. [2] it is smaller than rx- 
The blue line is a linear interpolation of such points. 



In the current literature (see e.g. Refs. [5l[l8l[Zl[ni[26l[29l[2Zl[T9l[20l[ai[l2l 
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I13[ I22p there are no analogous treatments, since it is always assumed an adiabatic 
dynamics instead of a fully radiative one and only the latest afterglow pha ses are 
addressed. It is usually assumed the Blandford-McKee self similar solution^QI for the 
adiabatic dynamics 7 cx r~^/^. A critical analysis of the applicability to GRBs of 
this approximate dynamics, instead of the exact solutions in Eqs. (l5l6l7p . has been 
presented in Ref. [1] as recalled in the introduction. The most widely applied formula 
for the EQTS apparent radius in the above mentioned current literature is the one 
proposed by SarM 

rx = 3.91 X IQis {E52/nif^ [Tdays/{1 + z)]'^"" cm, (12) 

where E52 is the initial energy of the shell in units of 10^^ ergs, ni is the CBM 
density in units of 1 particle/cm'^, T^ays is the arrival time at the detector of the 
radiation measured in days, and z is the source cosmological redshift. Waxman 
et al.J^ however derived a numerical factor of 3.66 x 10^^ instead of 3.91 x 10^^. 
Eq. ([T2l) cannot be compared directly with Eq. (fTT|) . since they assume two different 
dynamical regimes. Therefore, we first plot in Fig. |4] the exact solution for r±^ given 
by Eq. ([Tl]) both in the fully radiative case, using Eqs. (|2l3| ). and in the adiabatic 
one, using Eqs. ()5|6| ). We see that they almost coincide during the prompt emission 
while they start to diverge in the early afterglow phases, following the behavior of 
the corresponding Lorentz gamma factors (see Fig. [T] and details in Ref. Jj). 

Now, in Fig. [5] and for the same initial conditions assumed in previous figures, 
we plot the maximum transverse EQTS visible radius rj_ given by Eq. ljlip in the 
adiabatic case, using Eqs. (|5|6|) . together with Eq. ([T2|) . proposed in Ref. |5l and with 
the corresponding modification in the numerical factor proposed in Ref. 1181 From 
such a comparison, we can see that the approximate regime r±_ cx t^J^ overestimates 
the exact solution for during the prompt emission and the early afterglow phases. 
It is asymptotically reached only in the latest afterglow phases {ta > lO'^ s), in 
the very small region in which the approximate power-law dynamics starts to be 
applicable (see details in Ref. |T|). However, there is still a small discrepancy in 
the normalization: the constant numerical factor in front of Eq. p^ should be ^ 
3.1 x 10^^ instead of 3.91 x 10^^ or 3.66 x 10^^ to reproduce the behavior of the exact 
solution for large ta- Moreover, we must emphasize that, in analogy with what we 
obtained in the fully radiative case (see Fig. |3]), in the early phases of the prompt 
emission, when the fireshell Lorentz 7 factor is almost constant, the maximum 
transverse EQTS visible radius r± given by Ea. (|TT|) can only be considered an 
upper limit to the actual value of the EQTS apparent radius in the sky. In such 
phases the approximation implied by Eq. (jl2p can therefore be even worse that what 
represented in Fig. [S] 

5. Conclusions 

Within the fireshell model, using the exact analytic expressions for the fireshell equa- 
tions of motion and for the corresponding EQTSs in the fully radiative condition. 
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Arrival time (tg) (s) 

Fig. 4. Comparison between the maximum transverse EQTS visible radius rx computed with 
Eq. mil l in the adiabatic (blue line) and in the fully radiative (red line) cases respectively. 

we analyzed the temporal evolution of the distribution of the extended afterglow lu- 
minosity over the EQTS during the prompt emission and the early afterglow phases. 
We find that, at the beginning of the prompt emission [ta — 5 seconds), when 7 is 
approximately constant, the most luminous regions of the EQTS are the ones along 
the line of sight. As 7 starts to drop {ta ^ 30 seconds), the most luminous regions 
of the EQTSs become the ones closest to the boundary of the visible region. This 
transition in the emitting region may lead to specific observational signatures, i.e. 
an anomalous spectral evolution, in the rising part or at the peak of the prompt 
emission. The EQTSs of GRB extended afterglows should therefore appear in the 
sky as point-like sources at the beginning of the prompt emission but evolving after 
a few seconds into expanding luminous "rings" , with an apparent radius evolving 
in time and always equal to the maximum transverse EQTS visible radius r_L. 

We also derive an exact analytic expression for r_L, both in the fully radiative 
and in the adiabatic conditions, and we compared it with the approximate formulas 
commonly used in the current literature in the adiabatic case. We found that these 
last ones can not be applied in the prompt emission nor in the early afterglow 
phases. Even when the asymptotic regime is reached {ta ^ 10^ s), it is necessary a 
correction to the numerical factor in front of the expression given in Eq. ()12[) which 
should be - 3.1 X 10^^ instead of 3.91 x 10^^ or 3.66 x 10^^ 
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10° 10^ 10^ 10^ 10"* 

Arrival time (tg) (s) 



Fig. 5. The EQTS apparent radius in the sky as a function of the arrival time ta in the adiabatic 
regime. The red line represents the maximum transverse EQTS visible radius rx following Eq.JTT}. 
The green line represents Eq. lll2l l proposed by Sari. The blue line represents the corresponding 
modification in the numerical factor proposed by Waxman et al. These last two lines are almost 
coincident on the scale of this plot. 

References 

1. C. L. Bianco and R. Ruffini, J 633 (2005) L13. 

2. P. Meszaros, Reports of Progress m Physics 69 (2006) 2259. 

3. P. Couderc, Annales d'Astrophysique 2 (1939) 271. 

4. M. J. Rees, Nature 211 (1966) 468. 

5. R. Sari, ApJ 494 (1998) L49. 

6. A. Panaitescu and P. Meszaros, ApJ 493 (1998) L31. 

7. J. Granot, T. Piran and R. Sari, ApJ 513 (1999) 679. 

8. C. L. Bianco, R. Ruffini and S.-S. Xue, A&A 368 (2001) 377. 

9. C. L. Bianco and R. Ruffini, ApJ 605 (2004) LI. 

10. C. L. Bianco and R. Ruffini, ApJ 620 (2005) L23. 

11. A. Gruzinov and E. Waxman, ApJ 511 (1999) 852. 

12. Y. Oren, E. Nakar and T. Piran, MNRAS 353 (2004) L35. 

13. J. Granot, E. Ramirez-Ruiz and A. Loeb, ApJ 618 (2005) 413. 

14. Y. F. Huang, K. S. Cheng and T. T. Gao, ApJ 637 (2006) 873. 

15. Y.-F. Huang, Y. Lu, A. Y. L. Wong and K. S. Cheng, Chin. J. Astron. Astrophys. 
(ChJAA) 7 (2007) 397. 

16. E. Waxman, ApJ 491 (1997) L19. 

17. J. Granot, T. Piran and R. Sari, ApJ 527 (1999) 236. 

18. E. Waxman, S. R. Kulkarni and D. A. Frail, ApJ 497 (1998) 288. 



November 10, 2011 1:9 WSPC/INSTRUCTION FILE EQTS'Lum 



The luminosity evolution over the EQT3 in GRBs prompt emission 11 

19. T. J. Galama, D. A. Frail, R. Sari, E. Berger, G. B. Taylor and S. R. Kulkarni, ApJ 
585 (2003) 899. 

20. J. Granot and A. Locb, ApJ 593 (2003) L81. 

21. G. B. Taylor, D. A. Frail, E. Bcrgcr and S. R. Kulkarni, ApJ 609 (2004) LI. 

22. J. Granot, MNRAS 390 (2008) L46. 

23. D. A. Frail, S. R. Kulkarni, L. Nicastro, M. Feroci and G. B. Taylor, Nature 389 
(1997) 261. 

24. G. B. Taylor, E. Momjian, Y. Pihlstrom, T. Ghosh and C. Salter, ApJ 622 (2005) 
986. 

25. Y. M. Pihlstrom, G. B. Taylor, J. Granot and S. Doclcman, ApJ 664 (2007) 411. 

26. P. M. Garnavich, A. Loeb and K. Z. Stanek, ApJ 544 (2000) Lll. 

27. B. S. Gaudi, J. Granot and A. Loeb, ApJ 561 (2001) 178. 

28. K. loka and T. Nakamura, ApJ 561 (2001) 703. 

29. J. Granot and A. Loeb, ApJ 551 (2001) L63. 

30. R. D. Blandford and C. F. McKee, Physics of Fluids 19 (1976) 1130. 

31. R. Ruffini, C. L. Bianco, P. Chardonnet, F. Fraschetti and S.-S. Xue, ApJ 555 (2001) 
L113. 

32. R. Ruffini, A. G. Akscnov, M. G. Bcrnardini, C. L. Bianco, L. Caito, P. Chardonnet, 
M. G. Dainotti, G. de Barros, R. Guida, L. Izzo, B. Patricelli, L. J. R. Lemos, M. Ro- 
tondo, J. A. R. Hernandez, G. Vereshchagin and S. Xue, The Blackholic energy and 
the canonical Gamma^Ray Burst IV: the "long," "genuine short" and "fake-disguised 
short" GRBs, in XIII Brazilian School on Cosmology and Gravitation, eds. M. Novello 
and S. Perez Bcrgliaffa, American Institute of Physics Conference Series, Vol. 1132 
(2009), pp. 199-266. 

33. R. Ruffini, C. L. Bianco, P. Chardonnet, F. Fraschetti and S.-S. Xue, ApJ 581 (2002) 
L19. 

34. R. Ruffini, C. L. Bianco, P. Chardonnet, F. Fraschetti, V. Gurzadyan and S.-S. Xue, 
IJMPD 13 (2004) 843. 

35. T. Piran, Phys. Rep. 314 (1999) 575. 

36. C. L. Bianco and R. Ruffini, ApJ 644 (2006) L105. 



